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Taking into account effects of late energy injection, we examine big bang nucleosynthesis (BBN) 
constraints on axino dark matter scenarios with long-lived charged sleptons. We calculate 4-body 
slepton decays into the axino, a lepton, and a quark-antiquark pair since they govern late hadronic 
energy injection and associated BBN constraints. For supersymmetric hadronic axion models, we 
present the obtained hadronic BBN constraints and show that they can be more restrictive than the 
ones associated with catalyzed BBN via slepton-bound-state formation. From the BBN constraints 
on hadronic and electromagnetic energy release, we find new upper limits on the Peccei-Quinn scale. 



PACS numbers: 98.80.Cq, 95.35.+d, 12.60.Jv, 95.30.Cq 

I. INTRODUCTION 

It is widely believed that the Standard Model of par- 
ticle physics is not complete and that new physics is re- 
quired for a compelling theoretical description of nature. 
The Large Hadron Collider (LHC) may give us a first 
opportunity to test models of new physics directly. 

In addition to collider experiments, the early Universe 
offers many ways to probe new physics. In particular, the 
observed dark matter density and the framework of big 
bang nucleosynthesis (BBN) impose valuable restrictions 
on the parameters of new physics models. 

Supersymmetry (SUSY) and the Peccei-Quinn (PQ) 
symmetry are particularly well-motivated extensions of 
the Standard Model each of which comes with a com- 
pelling dark matter candidate. If both of these extensions 
are realized simultaneously, significant contributions to 
the cold dark matter density f2dm can reside in the axion 
and/or the lightest supersymmetric particle (LSP). 

An attractive realization of this scenario with con- 
served R-parity has the axino a — the fermionic part- 
ner of the axion — as the LSP and a charged slep- 
ton ^1 as the next-to-lightest supersymmetric particle 
(NLSP). In fact, the axino is an extremely weakly inter- 
acting particle and a promising dark matter candidate 
beyond the minimal supersymmetric Standard Model 
(MSSM) [ll,[l,i,i,[ES[7|,[i. Since the axino LSP allows 
for a long-lived charged slepton — such as the lighter stau 
Ti — which should be easy to discover at the LHC, those 
scenarios are appealing not only from a theoretical point 
of view but also from a phenomenological one. 

The requirement of avoiding overclosure of the Uni- 
verse from an overly efficient production of axino dark 
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matter allows one to derive upper limits on the reheat- 
ing ternperature [3, H, 0, S, M, [13, [HI- In our recent 
Letter [8| , we have shown how those constraints become 
more restrictive in the above scenarios once catalyzed 
BBN (CBBN) and the associated production of primor- 
dial lithium-6 and beryllium-9 [13, 0, [3] are consid- 
ered. Moreover, while axion physics provides us with a 
lower limit [H, M, [13 



/a > 6 X 10** GeV , 



(1) 



we found that CBBN also imposes upper limits on fa 
that depend on the slepton mass [8|. 

In this Letter we further develop the BBN constraints 
on axino LSP scenarios with long-lived charged sleptons 
by considering hadronic and electromagnetic energy in- 
jection. If the li is relatively long-lived, decay modes 
with hadrons in the final state can affect BBN and are 
thus strongly restricted. We analyze the most important 
contribution of this kind: the hadronic 4-body decay of 
li into the associated lepton, the axino, and a quark- 
antiquark pair. We also include the electromagnetic en- 
ergy injection, which is derived from the 2-body decay of 
the h into the the associated lepton and the axino. In 
fact, the BBN constraints on hadronic and electromag- 
netic energy release allow us to derive new upper limits on 
fa which appear in addition to those imposed by CBBN. 



II. PARTICLE PHYSICS SETTING 

We consider SUSY hadronic (or KSVZ ^\^) axion 
models .20:] in which the interaction of the axion multiplet 
<f> with the heavy KSVZ quark multiplets Qi and Q2 is 
described by the superpotential 



WpQ = y$0iQ2 



(2) 



with the Yukawa coupling y and the quantum numbers 
given in Table [D With the 2-component fields listed in 



2 



TABLE I: Quantum numbers of the axion multiplet $ and the 
heavy KSVZ quark multiplets Qi^2 considered in this work. 



chiral muhiplet 


U(1)PQ 


(SU(3)c,SU(2)l)y 




+1 


(l,l)o 




-1/2 


(3,l)+eQ 


Q2 = Q2 + V2q20 + F2ee 


-1/2 


(3*,l)-eQ 



that table, the axino and the heavy KSVZ quarks are 
described by the respective 4-coniponent fields, 



^ and g= M^^ 



(3) 



For the heavy KSVZ (s)quark masses, we use the SUSY 
limit Mg^ ^ = Mq = y{(f)) ~ yfa/V^ with both y and (0) 

taken to be real by field redefinitions.^ Accordingly, the 
phenomenological constraint (H)) implies, for y — 0(1), 
a large mass hierarchy between the KSVZ fields and the 
weak and the soft SUSY mass scales, 



III. ENERGY RELEASE IN STAU DECAYS 

Decays of the tr NLSP into the a LSP are associated 
with the emission of Standard Model particles. If in- 
jected at cosmic times t > 100 s, i.e., during or after 
BBN, those particles can reprocess the primordial light 
elements significantly via hadrodissociation and/or pho- 
todissociation. At earlier times, energetic particles are 
stopped efficiently through electromagnetic interactions 
so that hadrodissociation/photodissociation becomes an 
important issue only for > 100 s; see Ref. ,25j 
and references therein.'^ 

In this section, we explore the (average) electromag- 
netic/hadronic energy emitted in a single tr decay: 
ecm/had- Wc wiU usc this quantity below to investigate 
whether successful BBN predictions are preserved. 

The electromagnetic energy release eem is governed by 
the tau emitted in the main decay channel tr — » ra with 
an energy Er — (to~ — m|-|-m^)/(2m^:) in the rest frame 
of the tr . For stau decays at cosmic times t > 100 s, 
the emitted tau decays before interacting elect romagnet- 
ically. As each r decays into at least one v, which does 
not interact electromagnetically, only a fraction of Ej- 
contributes [28, ^31 . We use the conservative estimate 



(4) 



While our study can easily be generalized, we focus on 
the case in which the NLSP is a purely right-chiral lighter 
stau, Ti = Tr, which is a good approximation at least for 
small tan /3. Its coupling to ncutralinos is then dominated 
by the bino coupling. For further simplicity, we assume 
that mixing in the neutralino sector is such that one of 
the neutralino states is an (almost) pure bino. In fact, our 
cosmological considerations rely on a spectrum in which 
that state is the lightest neutralino, Xi = B, while our 
results for the tr decays are not restricted to this case. 
We treat the axino mass toj as a free parameter which is 
bounded from above by the stau NLSP mass rrir- 

In the considered R-parity conserving setting, the life- 
time of the Tr NLSP is governed by the decay tr — > to. 



Tr = 1/r: 



tot 



i/r(m 



(5) 



where the dominant contributions to the partial width 
r(TR — > to) occur at the two-loop level. ^ Using a heavy 
mass expansion in powers of 1/ fa [2l| . we have recently 
calculated the leading term that governs ^ [1]. 



To avoid ambiguities related to the fa definition, fa = \/2{<j>), the 
effective axion interactions obtained from ((2]| after integrating 
out the heavy KSVZ fields are given in Sect. 3 of Ref. g. 
As noted in H [1, HHH, the 3-body decay Tfi — > Ta7 occurs 
already at the 1-loop level, but for most of the phase space is 
subdominant. Moreover, compared to the 3-body decay, the 4- 
body decay modes Ta77, Tal~^l~, raqq are suppressed by 

an additional factor of a; cf. Fig. [T] below. 



0.3 Er = 0.3- 



2mr 



(6) 



to avoid that the electromagnetic BBN constraints pre- 
sented in Sect. IVII are overly restrictive. 

The hadronic energy release Chad is governed by the 
quark-antiquark pair emitted in the 4-body decay tr — > 
Taqq. Mesons from decays of the t's emitted in the main 
decay channe l tr — > ot typically decay before interacting 
hadronically [23[. In fact, since mesons typically decay 
before interacting with the background nuclei, only the 
nucleons originating from hadronization of the qq pair 
need to be taken into account for Tr > 100 s. We thus 
consider only qq pairs with an invariant mass TOqq above 



the mass of a pair of nucleons, 
when calculating 



2 GeV, 



Chad 



1 



- tot "'"'qH 



rfr(fR^T5qq) 



According to the value of Wqq, all quark flavors that can 
occur in the final state are taken into account in the 
differential decay rate dVij^ raqq) / dniqq. 

The Feynman diagrams of the dominant contributions 
to the 4-body decay tr Toqq are illustrated in Fig. [TJ 
where the blobs represent heavy KSVZ (s)quark loops; 
cf. Fig. 2 of Ref. [1] . As in the calculation of the 2-body 



^ The presence of additional slow hadrons can also affect BBN 
through proton-neutron interconversion processes |26| . The as- 
sociated limits are not considered since they are typically milder; 
cf. [23I for a discussion in G LSP scenarios with a stau NLSP. 
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FIG. 1: Feynman diagrams of the dominant contributions 
to the stau NLSP decay tr, ^ raqq in a SUSY hadronic 
axion modeL Here the hghtest neutrahno is assumed to be a 
pure bino Xi ~ B and the tan mass is neglected. The blob 
represents the loops involving the heavy KSVZ (s)quark fields 
with quantum numbers given in Table HI cf. Fig. 2 in Ref. Q]. 



decay [8|, we work in the limit nir which is justified 
since <^ rrir. Taking advantage of (|4]), we use again 
the method of heavy mass expansion to expand in heavy 
(s)quark masses, keeping only the leading terms oc 1/ fa 
in the calculation of the decay amplitude. 

Although the decay first occurs at the 1-loop level, 
the shown 2-loop diagrams are also important as their 
contribution to the amplitude is enhanced by large log- 
arithms ln{yfa/V^mr) ~ 20, e.g., for my/y = 100 GeV 
and fa = GeV. This enhancement largely compen- 
sates for the loop suppression factor governed by the fine- 
structure constant a = e^/ (47r). Therefore, only the lead- 
ing logarithmic term of the 2-loop contribution is kept. 

There are more diagrams at the 2-loop level than those 
shown as part of Fig. [T] One group has a 'y, Z line that 
splits into the qq pair after originating either from the 
Tr, or T in the loop or from one of the heavy (s)quarks 
in the blob. However contributions of these diagrams are 
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FIG. 2: Energy spectrum of the quark-antiquark pair emit- 
ted in the four-body decay tr — > raqq with an invariant mass 
rriqq, normalized to the total decay rate « r(TR —> to). 
The shown quantity is independent of eg. The other param- 
eters are set to fa = 10^'^ GeV and y = 1. 



canceling or not enhanced by large logarithms and there- 
fore neglected in our calculation. Similarly, there are 
diagrams with gluons (more than one because of color 
conservation) emitted by the heavy (s)quarks in the blob 
at higher order in the strong coupling, turning into nucle- 
ons. These are also highly suppressed and thus neglected. 



^ tot 



our results for the differential decay 
Taqqj/drriqq normalized to the total 
are shown for (ms, mr, 'nig) com- 



In Fig-^ H 
rate c?r(TR - 
decay rate 

binations of (10 GeV, 100 GeV, 1.1 to?), (50 GeV, 
150 GeV, 1.1 TO?), (10 GeV, 150 GeV, 1.1 to?), and 
(10 GeV, 150 GeV, 1.02 m?) by the solid, dash-dotted, 
dotted, and dashed curves, respectively. The shown 
quantity is independent of eg. The other parameters 
are set to fa — 10^^ GeV and y = 1. The kinks at TOqq ~ 
3 GeV and 10 GeV mark the respective thresholds for 
the production of cc and bb pairs. Moreover, one can see 
the contribution of the Z-boson resonance at mqq = Mz , 
which is present for Am = m^ji — 'ma — mr > Mz- As part 
of the integrand in ([7|) , the shown energy spectrum of the 
emitted qq pair is the crucial quantity in the calculation 
of Chad- Based on the numerical results from the integra- 
tion in ([7|) , we obtain the hadronic BBN constraints pre- 
sented in Sects. |V] and IVII We note that a more precise 
calculation of the constraints would require, in addition 
to dT(Tji — > raqq)/dmqq, a treatment of the fragmenta- 
tion of the quarks into hadrons and of the propagation of 
the resulting hadron spectra when computing the abun- 
dances of primordial light elements; cf. [23l. [25l. [30l. [3]| . 
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IV. COSMOLOGICAL SETTING 

In addition to and ecm/had: the stau NLSP yield 
prior to decay, = n:f^/s, is another quantity that is 
crucial for our study of cosmological constraints. Here s 
denotes the entropy density and n^^ = n~+ +n~- = 2 n~- 
the total TR number density prior to decay for an equal 
number of positively and negatively charged trS. In con- 
trast to and eom/had, Yr depends on the thermal his- 
tory of the early Universe. In this Letter we assume 
a standard thermal history with a reheating tempera- 
ture after inflation in the range: T^^^ < < fa, where 
^dec ~ 1^^/25 is the temperature at which the tr NLSP 
species decouples from the primordial plasma. This has 
a number of implications: 

(i) Focussing on a standard thermal history, we assume 
that effects of the saxion — the bosonic partner of 
the axino that appears in addition to the axion — 
are negligible. For a non-standard thermal history 
associated with significant late entropy production 
in saxion decays [H, [H, [H, HH] , cosmological con- 
straints [10] including those considered in this work 
can be affected depending on the saxion properties. 

(ii) For Tr < fa, no PQ symmetry restoration takes 
place after inflation. In fact, we assume that the 
PQ symmetry was broken before inflation and not 
restored afterwards. For large fa such that axions 
are never in thermal equilibrium with the primor- 
dial plasma, the relic axion density fla is then gov- 
erned by the initial misalignment angle 8^ of the 
axion field with respect to the CP-conserving posi- 
tion; cf . [l^, 113] and references therein. This allows 
us to keep the presented constraints conservative by 
assuming fla ^ ^dm which is possible even for fa 
above 10^'* GeV since 9^ can be sufficiently small. 

(ni) For TJJ^ < Tr, the tr NLSP decouples as a weakly 
interacting massive particle before its decay into 
the axino LSP so that F^r is the thermal relic stau 
abundance, which does not depend on Tr. 

The thermal relic stau abundance prior to decay can 
be calculated numerically. Its value depends on details 
of the SUSY model such as the mass splitting among the 
lightest Standard Model superpartners [36| or the left- 
right mixing of the stau NLSP [13, [1^. With the focus 
on the Tr NLSP setting in this Letter, we consider three 
characteristic approximations: 

~ «: X 10-12 (^^) , = 0.7, 1.4, 2.8, (8) 

where n accounts for typical differences in the annihila- 
tion processes of the sleptons. The value k = 0.7 corre- 
sponds to the case with rrig = 1.1 and rrir «C m^^ji, in 
which primordial stau annihilation involves only staus in 



the initial state [s^.'' The yield associated with k = 1.4 
is encountered if there is either additional stau-slepton 
coannihilation corresponding to < m^.ji < 1.1 (36j 
or additional stau-bino coannihilation corresponding to 
TOy < mjj < 1.1 (cf. contours close to the dashed 
line in the right panel of Fig. 3 in Ref. 39]). For an 
approximate degeneracy of rrir with both me^ji and mj^, 
simultaneous stau-slepton-bino coannihilation can lead 
to an even larger k = 2.8 in ([8]) (cf. contours close to 
the dashed line in the left panel of Fig. 3 of Ref. [39]) . 

A non-thermally produced (NTP) axino density[ll, d, 
ii 

nr^'h^ ^maYrs{To)h^/pc (9) 

emerges since one axino LSP is emitted in each stau 
NLSP decay; pc/[s{To)h^] = 3.6 x lO-^GeV [H]. This 
density must not exceed the dark matter density [13] 

f^I-/i2 = o.io5l°:02i^ (^0) 

where a nominal 3cr range is indicated and where 
h = 0.73lg;o3 denotes the Hubble constant in units of 
100 kmMpc~is~i. Additionally taking into account the 
thermally produced axino density flT^ and the ax- 

ion density fla, one obtains the dark matter constraint 
on the stau abundance prior to decay: < Y^Xm '^ith 

.„-..(2i^)(ii^).,. 

For the most conservative case + ila ^ ^dm, we 
have illustrated the constraint on ms and rn? obtained 
by confronting Y^ with pT|) already in Ref. [g]. 
For T^r > 10'^ s, additional upper limits on Y— — 

occur since the negatively charged t^s can form (^He ) 
and (^BcT^) bound states and can thereby catalyze pri- 
mordial ^Li and ^Be production in excess of observation- 
ally inferred limits [l3,|il[ll.'^ In fact, the CBBN con- 
straints on rria and m^f obtained by confronting Y^ with 
the CBBN-induced T^-dependent upper limits given in 
Ref. [13] have been presented for the first time in Ref. [^ . 

V. HADRONIC BBN CONSTRAINTS 

Hadronic energy injection can affect the abundance 
of primordial deuterium substantially via hadrodisso- 
ciation of helium-4. In fact, we focus on the con- 
straint on hadronic energy release imposed by the pri- 
mordial abundance of D in this section. While addi- 
tional constraints on late energy injection are imposed 



^ The bino mass nij^ = 1.1 considered in Ref. l36l represents a 
typical mass splitting in regions with rUg > m^: encountered in 
scenarios such as the constrained MSSM (CMSSM). 

^ We note that Ref. [4lj has questioned the efficiency of the cat- 
alyzed production of ^Be obtained in Refs. [Tsl . Il4j . and an- 
nounced a further clarification of this point in the future. 



5 



> 



0.5 



0.2 



0.1- 



(a) 



10^ 


s 




5 














-V— \ 






■ V- 




2 


5 X 


10^ s 








T-)SCV . ■ / 


S 0.5 


10'^ 


s 






fa = 

\eQ\ = 


3 X 1013 GeV 
1/3 


/i 

/J" 


0.2 


rrig = 
Yf = 


1 . 1 rrif 

0.7 X 10-12 (mf/lTeV) 




0.1 


1 


10 


100 1000 




ma [GeV] 




(b) 



!a = 3 X 10^2 GeV 
leg! = 1/3 



Yf = 2.8 X 10-12 (mf /lTeV)\ 



10^ s 



5 X 10^ s 




10'"' s 



10 100 
ma [GeV] 



1000 



FIG. 3: Cosmological constraints on the masses of the a LSP and the tr NLSP for (a) /a = 3 x IQi'^ GeV, = 1.1 m^, 
Yr given by (HJ with k = 0.7 and (b) = 3 x lO^^ GeV, = 1.02 m^=, Yf given by (O with k = 2.8. In both panels, 
leg I = 1/3 and y = \. The hadronic BBN constraints associated with (I13p and (|14p disfavor the regions enclosed by the 
respective short-dash-dotted (blue) lines. The CBBN constraints associated with (|15|l and (|16|) disfavor the regions below the 



long-dash-dotted (red) and the solid lines. Contours of = 10 , 5 x 10 , and 10 s are shown by the dotted lines. On (above) 
the gray band, i^^^'' G fldm (S^z^^^ > 0.126). The region with mg > mf is not considered as we focus on the 5 LSP case. 



by the primordial abundances of ''He, '^He/D, ''Li, and 
^Li/'^Li [23, 24, 25, 42,, 44, 4^, the hadronic D con- 
straint is the dominant one in the region allowed by the 
CBBN constraints, Tt, < 10'^ s. This can be seen, e.g., in 
Figs. 38-41 of ReL [23| and in Figs. 6-8 of Ref. [2|. Also 
the D constraint on electromagnetic energy injection is 
not considered in this section since it becomes relevant 
only for Tf^ > IC* s which is already disfavored by CBBN; 
see Sect. IVII for details. 

To derive the hadronic BBN constraint imposed by D, 
we use the TNLSP-dependent upper limits (95% CL) on 

6iad = Chad i^NLSP (12) 

given in Fig. 9 of Ref. as obtained in Ref. [23] for 
observationally inferred primordial D abundances of 

D/H|„,ean = (2.78l°;^^) X 10-5 (severe) , (13) 
D/H|high = (3.98t^j?) X 10-5 (conservative) . (14) 

Here tnlsp = t?, Inlsp = i?, and ehad is given in Q. 
With the calculated ehad 7 upper limits on can be 
derived from the considered upper limits on ^had [4^ : 
^?Tad = eSVehad. By confronting Y^ with Fi^, we 
then obtain the regions in the parameter space that are 
disfavored by the hadronic BBN constraints. 

In Fig. [3] the obtained hadronic BBN constraints are 
shown by short-dash-dotted (blue) Hues. The labels D^'g^ 
and D^"^*' indicate the constraints associated with (fT5)) 
and respectively. The regions enclosed by the cor- 
responding lines are disfavored by an excess of D above 



the respective observationally inferred abundance. The 
absence of a 0^°^*^ line in panel (a) and the difference be- 
tween the D^'g^ and DJ^° ^'^ lines in panel (b) indicate the 
sensitivity on the observationally inferred D abundance. 

We also show the CBBN constraints associated with 
the observationally inferred limits on the respective pri- 
mordial fractions of ^Li g^,!!!,!!!] and ^Be [11, 

"^Li/Hlob. < 10-"-10-i° , (15) 
^Be/Hlobs < 2.1 X 10-13 ^ (16) 

as obtained by confronting Y~- with the limits in Fig. 5 

of Ref. [ll . The range is indicated by pairs of long- 
dash-dotted (^Li, red) lines and (|16p by solid (^Be) lines. 
The regions below those lines are disfavored by an excess 
of ^Li and ^Be above the respective limits. The region 
with f^j"'"^ S indicated by the gray band, where 

the region above that band is the one excluded by (fTT|) in 
the conservative case with +ria ^ f^dm- The dotted 
lines are contours of = lO'^, 5 x lO'^, and 10^ s. 

For [eg I = 1/3 and i/ = 1, we consider /a = 3 x 
10" GeV, = l.lrrir, Yr given by ^ with n = 0.7 
in panel (a) and fa = 'i'x 10^^ GeV, nig — 1.02 m?^, V? 
given by ([8]) with k — 2.8 in panel (b). In both panels, 
one finds my values disfavored by the hadronic BBN con- 
straints that are significantly larger than the ones disfa- 
vored by the CBBN constraints. However, for jegl — 1/3, 
y = 1, rrig = 1.1 my and k — 0.7, even the more re- 
strictive 'severe' hadronic BBN constraint disappears for 
fa ^ lO"'^ GeV. In contrast, the CBBN constraints re- 
main until fa < lO'^^ GeV — cL Fig. [D^a) — given the limit 
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(a) ma[GeV] (b) ma[GeV] 

FIG. 4: Cosmological constraints on the masses of the o LSP and the tr NLSP for (a) rrig = 1.1 m^^ and (b) rrig = 1.02 m^=. 
In both panels, fa = 10^^ GeV, \eQ\ = 1/3, y = 1, and is given by (O with k — 1.4. The BBN constraints, the dark matter 
constraint, the contours, and the other regions are indicated as in Fig. [31 



> 80 GeV [T^ from searches for long-lived staus at 
the Large Electron Positron (LEP) collider. 

For rUg = 1.02 m^r and enhanced by stau-slepton- 
bino coannihilation such that k = 2.8, Fig. [3jb) shows 
that 'severe' and 'conservative' hadronic BBN constraints 
are still encountered at /a = 3 x 10^^ GeV. This demon- 
strates the sensitivity of those constraints on the stau 
NLSP yield. In fact, the restrictive I^^^^ constraint re- 
mains until fa < 10^2 GeV for k = 2.8; cf. Fig.IHb). 

In the latter case, more restrictive hadronic BBN con- 
straints are encountered not only due to the enhanced 
yield but also due to the smaller bino mass which is asso- 
ciated with larger values of ehad and thereby with more 
restrictive Y~^^^ values. Fig. d] shows this more clearly 
for fa = 10^3 GeV, leg] = 1/3, y = 1, and k = 1.4, where 
constraints and values are indicated as in Fig. [31 The 
two panels seen here only differ in the bino-stau mass 
ratio. The one on the left illustrates the case where 
we have — 1.1 m^: and stau-slepton coannihilation 
while the panel on the right illustrates the case where 
we have rUg = 1.02 m^; and stau-bino coannihilation, 
where the different origins of coannihilation only serve to 
give the same in the two cases. However, this change 
in bino-stau mass ratio alters the result sufficiently to 
produce much more restrictive D constraints in the lat- 
ter case. We can trace this effect back to ehad through 
Fig. [21 which illustrates the considerably higher curve for 
the case of rrig — 1.02 as compared to = l.lm^, 
with the same ms, m?, and fa values (the top two curves 
in the plot). A thorough study of this and of other as- 
pects of the hadronic BBN constraints will be presented 
in a forthcoming publication (46j . Our objective in the 
present Letter is (i) a first presentation of the hadronic 



BBN constraints in the case with the axino LSP and a 
charged slepton NLSP and (ii) to demonstrate explic- 
itly that they can be more restrictive than the associated 
CBBN constraints. 

Before proceeding let us comment on the potential in- 
terplay between late energy injection and CBBN. The 
CBBN limits adopted from Ref. [T^ have been derived 
for an abundance of D obtained with standard BBN. 
For an increased D abundance from hadrodissociation 
of "^He, CBBN of ^Li and ^Be becomes more efficient. 
This is evident for ^Li since its catalysis proceeds via 
(^HeTp ) -I- D — > ^Li -f Tj^ [l2| and since the primordial 
abundance of D stays significantly below the one of ^He 
at the relevant times (even for a maximum of observa- 
tionally tolerable hadrodissociation of ^He); cf. Fig. 2.4 
in Ref. An increased output of ^Bc results from the 
final step of its catalysis, {^BeT^)+n —^'^ Be-|-fj^ [T3l.lT3|. 
which becomes more efficient since an enhanced abun- 
dance of D increases the number of neutrons n at the 
relevant times [Slj; cf. Fig. 4 in Ref. [l^. Moreover, the 
debris of hadrodissociated ''He can hit ambient *He and 
thereby fuse additional ^Li [H [13, [3. The interplay 
of late energy injection and CBBN will thus lead to con- 
straints that can only be stronger than the ones presented 
in this Letter. Aiming at conservative limits, this allows 
us to neglect those intricacies which will have to be faced 
in future refinements of the presented constraints. 



VI. BBN CONSTRAINTS ON THE PQ SCALE 

In this section we show that the BBN constraints asso- 
ciated with hadronic and electromagnetic energy release 
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— l.lm^, Yr given by (|5} with k = 0.7 and (b) m^j = 1.02 m^, 
Yyf given by ((8)1 with k = 2.8. In both panels, m^/m~ <^ 1, \eQ\ = 1/3, and y = 1. The hadronic BBN constraints associated 
with (|13p and (|14p disfavor the regions in the upper right-hand corner enclosed by the respective short-dash-dotted (blue) lines. 
Electromagnetic BBN constraints associated with D disfavor the upper regions enclosed by the respective dashed (blue) lines 
and the ones associated with ^He/D the region above the double-dash-dotted (green) line. The regions above the long-dash- 
dotted (red) and the sohd hues are disfavored by the CBBN constraints associated with (|15|) and (|16|) . Contours of Tf — 10^, 
10'', and 10^ s are shown by the dotted lines. 



impose new upper limits on the PQ scale fa- In agree- 
ment with the results of the previous section, we find that 
those limits can be substantially more restrictive than the 
ones imposed by the CBBN constraints Q. 

Fig. [5] presents the new fa limits together with the 
CBBN limits for m|/m| < 1, jegj = 1/3, and y = 1. 
In panel (a) we have the generic case of rrig = l.lriir 
and Yr given by ^ with k = 0.7 while in panel (b) 
we have rrig = 1.02 m? and k — 2.8 representing the 
case with simultaneous stau-slepton-bino coannihilation. 
Contours of = 10^ 10^, and lO^s are shown by the 
dotted lines. Above the long-dash-dotted (red) and the 
solid lines, CBBN of ^Li and ^Be is in excess of the re- 
spective limits and p6p . The new hadronic BBN 
constraints associated with (|13p and (|14p disfavor the re- 
gions in the upper right-hand corner enclosed by the re- 
spective short-dash-dotted (blue) lines. Note that those 
constraints are provided only for t? > 100 s since we 
have not considered the typically milder limits associ- 
ated with proton-neutron interconversion processes p6j 
which become relevant for smaller [H, 0, . Nev- 
ertheless, the hadronic BBN constraints place limits on 
the PQ scale fa that become clearly more restrictive than 
the CBBN-induced limits towards large to? and/or large 
Yr- In fact, the hadronic BBN constraint on fa can be 
the dominant one already in a mass range, to? < 1 TeV, 
that is promising for a discovery of a long-lived stau at 
the LHC. 

While the above sets of BBN constraints correspond to 



the ones shown in the previous section (cf. Fig. [3]), we also 
indicate in Fig. [5] the electromagnetic BBN constraints 
imposed by primordial D and '^He/D. Our derivation 
of the electromagnetic BBN constraints proceeds as out- 
lined for the hadronic ones in Sect. |V]but relies on the 
conservative eom © Si,iid on upper limits on ^om = eom V?. 
Accordingly, we obtain the shown D^°^ and ■^He/D con- 
straints from the respective limits given in Fig. 42 of 
Ref. [2^ and the Dcm^ constraint from the respective 
limit given in Fig. 6 of Ref. [47|. Those Dcm^™"' and 
■^He/D constraints disfavor the upper regions enclosed by 
the respective dashed (blue) lines and the regions above 
the double-dash-dotted (green) lines in Fig. [5j 

Fig.Oshows that the electromagnetic BBN constraints 
appear only for t?^ > 10* s thereby excluding regions 
already disfavored by CBBN. Nevertheless, they support 
the finding that, e.g., values of the PQ scale at the scale 
of grand unification, fa ~ 10^^ GeV, will be in conflict 
with successful BBN in the considered scenarios once a 
long-lived charged slepton is observed at the LHC. 

Before closing let us discuss the robustness of the 
shown fa limits and address important sensitivities: 

• By considering ml/m- <C 1, the CBBN-imposed 
fa limits are conservative limits. Those constraints 
become more restrictive for ms — > to?. This is dif- 
ferent for constraints associated with late energy in- 
jection, where any bound can be evaded for a finely 
tuned m^-TO? degeneracy leading to Chad/om 0. 
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• The fa limits are sensitive to Y^. In settings with 
a sizable left-right stau mixing, an exceptionally 
small Yf is possible such that even the CBBN con- 
straints may be respected [13, HI] . 

• The fa limits depend on the quantum numbers of 
the heavy KSVZ fields. While ehad/om are indepen- 
dent of eg, T^: cx l/eg. The fa limits can thus be 
relaxed, e.g., by one order of magnitude for eg = 1. 

• The CBBN and hadronic BBN constraints in the 
case of the cr or piR NLSP are identical to the 
ones shown. The electromagnetic BBN constraints 
however will be more restrictive in the e^i NLSP 
case since all of the electron energy released in 
the Br NLSP decay will contribute: ecm = Eg. 

VII. CONCLUSION 

For axino LSP scenarios with a long-lived charged slep- 
ton NLSP, we have studied BBN constraints associated 



with hadronic and electromagnetic energy release. While 
the region with fa ^ 10"'^^ GeV is typically not affected, 
those constraints become significant for larger fa such 
that models with fa towards the grand unification scale 
are disfavored. The new BBN constraints on fa can be 
more restrictive than the recently obtained CBBN con- 
straints Q. This further tightens the upper limits on 
the reheating temperature discussed in Ref. [8] which are 
relevant for models of infiation and baryogenesis. 
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